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ABSTRACT 

In the context of brane- world scenarios, we study the effects produced by 
the exchange of virtual massive branons. A one-loop calculation is performed 
which generates higher-dimensional operators involving SM fields suppressed 
by powers of the brane tension scale. We discuss constraints on this scenario 
from colliders such as HERA, LEP and Tevatron and prospects for future 
detections at LHC or ILC. The most interesting phenomenology comes from 
new four-particles vertices induced by branon radiative corrections, mainly 
from four fermion interactions. The presence of flexible branes modifies also 
the muon anomalous magnetic moment and the electroweak precision ob- 
servables. 
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1 Introduction 



In recent years, it has been shown that a generic property of brane-world 
models [1] with low tension (r = f 4 <C A 4 , where r is the brane tension 
and A is the scale below which the description given by the brane-world 
scenario is appropriate) is the presence of new modes 7i- a (x) called branons 
which roughly correspond to excitations of the brane position along the extra 
compactified dimensions. The relevant tree-level phenomenology of branons 
has been studied for colliders and also for astrophysics and cosmology in 
terms of their mass M and brane tension parameter /, and it has been 
suggested that massive branons could be natural candidates for dark matter 
in this kind of models [2]. 

In this work we study the phenomenology of branon radiative corrections. 
Branon loops are interesting mainly for two reasons. First because preci- 
sion tests of the Standard Model (SM) usually enforce strong constraints on 
physics beyond it and thus make possible to reject many new models, or at 
least to set bounds on their parameters. The second reason is that branon 
loops provide new physical effects, such as four fermion interactions, which 
can be searched for in present and next generation colliders. 

As it is the case for branon tree-level effects, the loop corrections can be 
obtained from the effective action for branons described in detail in [3] . This 
effective action can be expanded in powers of dn/ f 2 and M 2 7r 2 // 4 [4, 5]: 



The zeroth order term is just a constant, the 0({n/ f) 2 ) contribution contains 
the branon free action: 



where M 2 ^ is the squared branon mass matrix corresponding to the different 
branon excitations ir a , with a running from one to the number of effective 
extra dimensions N. The couplings to the SM fields $ living on the brane 
(or any suitable extension of it) in the presence of a gravitational background 
which for simplicity we will assume to be flat (g^ = rj^), can be described 
at low energies by the action: 



SeffW = s%>[*] + s$M + s$M + ... 



(i) 






i 



+ ^(4^vr^-M> Q 7rV)T^ + 0(ir% (3) 



where T^ V M is the conserved SM energy-momentum evaluated in the back- 
ground metric: 



It is interesting to note that there is no single branon interactions due to 
the parity conservation on the brane by the gravitational action. Thus bra- 
nons are absolutely stable. This fact is crucial for the branon phenomenology 
and in particular for cosmology since it makes them natural WIMP candi- 
dates for dark matter [2] . In addition, the quadratic expression in (3) is valid 
for any internal extra-dimension space Kn, regardless the particular form of 
its metric j a ^. Indeed the form of the couplings only depends on the num- 
ber N of branon fields and the brane tension. Dependence on the geometry 
of the extra dimensions will appear only at higher orders. Here we are as- 
suming that the bulk D dimensional space-time (D = 4 + N) can be split 
as M. D = M 4 x K N , where M 4 is the standard four-dimensional Minkowski 
space and Kn is some compact and homogeneous space of dimension N with 
gaussian coordinates y a . Then the N branon fields (a = 1, ...,7V) can be 
chosen so that n a (x) = f 2 y a (x) where y a = y a (x) represents the position of 
the excited brane in the extra-dimension space K^. The brane ground state 
corresponds to n a = 0. 

From the action above it is clear that branons always interact by pairs 
with the SM matter fields. In addition, due to their geometric origin, those 
interactions are very similar to the gravitational ones since the 7r a fields 
couple to all the matter fields through the energy-momentum tensor and 
with the same strength suppressed by a f A factor. The interaction between 
bulk gravitons and SM fields is given by: 



where Mp = Mp/Air is the squared reduced Planck mass (Mp = 1.2 x 
10 19 GeV) and hffl are the Kaluza-Klein (KK) modes of the bulk graviton 
h^ u (x,y) corresponding to the 4x4 part of the bulk metric: 







9^u(x,y) = rj„ u + 




(6) 
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where M^~ 2 = Mp~ 2 /(Air) and Md is the D dimensional Planck scale (fun- 
damental scale of gravity) related with the usual four-dimensional Planck 
scale by Mj> = V(K N )M£~ 2 with V(K N ) being the volume of the internal 
space K N (notice that M 4 = M P ). By introducing a complete orthonormal 
set of functions f p = f p (y) on K N with normalization: 



and /p(0) = 1, the KK mode decomposition for the graviton field becomes 



When computing radiative corrections, divergent integrals appear. As 
our effective actions are not renormalizable all our results will be given in 
terms of some energy cut-off A, which could be taken as the value where the 
whole brane-world picture breaks down and a more fundamental approach is 
needed. Then our results will be given in terms of four parameters, namely 
the number of branons or extra dimensions N, the branon mass M (for 
simplicity we will assume at the end that all of them are degenerate M a p = 
5 a pMp = SapM), the brane tension scale / (r = f A ) and the cut-off A. 

The plan of the paper goes as follows: In section 2 we reobtain the result 
of [7] concerning the suppression of the coupling between SM fields on the 
brane and bulk fields, by integrating out the branon fields instead of using 
arguments based on normal ordering. In section 3 we study the effects of 
branon loops on the SM particle parameters and find the effective action 
describing the new induced interactions. The corresponding phenomenolog- 
ical consequences are considered in section 4, where we also set the bounds 
coming from the branon loops on the parameters /, M, iV and the scale A. 
Further constraints can be obtained from two loops effects and their impact 
on the electroweak precision observables and the muon anomalous magnetic 
momentum which can be found in section 5. In section 6 we summarize and 
comment our results and in Appendix A, B and C we define the divergent 
integrals appearing in our computations, the Feynman rules corresponding 
to the effective Lagrangian describing the branon loops effects, and the as- 
sociated cross-sections. 
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(8) 
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2 Graviton coupling suppression 



Probably the most immediate effect of virtual branons is the suppression of 
the coupling of SM particles and the KK modes bulk fields like the graviton. 
When branon fluctuations are taken into account this effective coupling is 
described by the action: 




due to the fact that the brane is no more sitting at tt = but is moving 
around this point. Now the branons fields can be integrated out in the usual 
way to find: 

Sh = t^E / <Pzh$(x)T? M (x)(fM) (10) 

where the f p expectation value is given by 

(f p (n))= f[dn]e^f p (n) (11) 

In the limit of massless branons, the branon effective action is just a non 
linear sigma model (NLSM) based on a coset which is isomorphic to K N . 
Therefore the invariant path integral measure should include an additional 
factor proportional to the square root determinant of the coset metric to 
ensure that quantum corrections do not spoil the Ward identities of the 
NLSM. The extra term in the measure amounts to an extra term in the 
effective action proportional to A 4 . This term is important when dealing 
with branon loop corrections to the branon self-interactions (for instance 
branon-branon elastic scattering) [6]. However in this work we are mainly 
interested in interactions between a couple of branons and SM particles and 
hence we can safely neglect this measure term. 

To compute the path integral above, we need to know the precise form of 
the f(n) functions which depends on the K N geometry. For example for the 
case of the torus K N = T N : 

fn(y) = exp (i^j (12) 

where n = (rii, ri2, ...njv) is a N dimensional vector with integer and positive 
or zero components and R is the torus radius (common for all coordinates). 
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Then we have 



(/«(*■)) = (exp 



nir 



RP 



N 



) = exp 



2R 2 f 



E»« G «*(°) 



(13) 



where G a p{x) is the branon propagator 

G«/j(:r-3/) = j dqe'^^ 



q 2 - M 2 + ie 



(14) 



and c/g = d 4 q/(2Tr) i . Using a cut-off A to regularize the divergent integral 
we find 



16tt 2 

Then the effective action becomes 

1 



A 2 - M 2 a log 



A 2 

M 



s h 



Mr 



E / ^xgnh^{x)T^ M {x) 



(15) 



(16) 



In other words, the effect of branon quantum fluctuations amounts to intro- 
ducing the KK mode dependent couplings which are given for toroidal 
compactification by 



A 2 



9n = exp 



327T 2 /? 2 / 4 i 



N 

E 



M 2 
l"#log 



' A^ 

M 



(17) 



Thus the coupling of SM matter to higher KK modes is exponentially sup- 
pressed. This result was first obtained in [7] for massless branons by using an 
argument based on normal ordering. Our derivation here is more natural in 
the context of the path integral treatment of branon quantum fluctuations, 
and also it can be applied to massive branons in any extra dimension space 
K N . In any case this coupling suppression has very interesting consequences 
from the phenomenological point of view. It improves the unitarity behavior 
of the cross section for producing gravitons from SM particles and, in ad- 
dition, it solves the problem of the divergences appearing even at the tree 
level when one considers the KK graviton tower propagators for dimension 
equal or larger than two. Moreover whenever we have v = Rf 2 A, KK 
gravitons decouple from the SM particles, so that at low energies the only 
brane-world related particles that must be taken into account are branons. 
In the following we will assume this to be the case and accordingly we will 
deal only with SM particles and branons. 
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3 Branon-loops effects on SM particles 

In order to study the effect of virtual branons on the SM particles, it is useful 
to introduce the SM effective action r^[$] obtained after integrating out 
the branon fields: 

e <>] = J [ dn ] e iS S Ml*] = e '/^sM (Det [ ])-i/2 (18) 

where the O operator is defined as O = A + B with: 

A a p(x,y) = -6 aP [d^ + M 2 a }5(x - y). (19) 



B aP {x-y) = jfSap 



M 2 



6(x-y), 



(20) 



This SM effective action can be computed in a systematic way by using 
standard procedures (see for instance [8]). Thus we have: 



and 



r e /^] = fd'xCsM + ^Tr(HO}), 



Tr(ln[0]) = Tr(ln[A]) + Tr(ln[l + BA ]). 



(21) 



(22) 



The first term does not depend on the SM fields and it can only contribute 
to the renormalization of the the cosmological constant. Expanding the log- 
arithm, we obtain the usual expression: 



r eff 
1 SM 



[$] = Jd 4 x£ 



z k=i K 



-Tr(BA 



-l\k 



/oo 
A£ SM + E r(t) [ $ ]. 
fc=l 



(23) 



(24) 



where A is the branon (n a ) free propagator: 

A Zp( x >y) =G a p(x-y). 
Then the first contribution reads: 

= l -jd i xd i yB^(x,y)G aP {x-y) = C 1 jd i xT s » M ^ (25) 
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where, assuming all the branons to be degenerate (M a = M, a = 1, . . . , N), 
the C\ constant is given by: 

= q *-M> + ie - (26) 

The second contribution to the effective action is: 

r< 2 >[$] = ~ t -ld 4 xd 4 yd 4 zdHB a(3 (x,y)G^(y~z)B y5 (z,t)G 5a (t-x) 
N r 

= -i— / d 4 xd 4 yd 4 pe- tp(y ' x) T' IU (x)T pa (y) (27) 
M 2 M 4 

where from now on T pv = Tglj and the integrals J {1) are defined in Appendix 
A. It is convenient to split the final expression into a local divergent term 
T^ 2) [$] and a non-local finite term: r^ L [$]. 

= ifpl + rgip]. (28) 

Now by using the equation of motion at the zero order: d^T^ v M = 0, the local 
piece can be written in terms of six constants Wi, (i = 1, 2, ..6): 

if } [$] = Jdx {W^T^ + w 2 thx 

+ W 5 T^U 2 T^ + W^U 2 T V V } (29) 
and the non-local one in terms of two functions Dj(p), (j = 1,2): 

= J dx dy dpe- ipiy - x) 

{D 1 {p)T^(x)T, l/ (y) + D 2 (p)T»(x)T»(y)}. (30) 

These two functions are given by: 

DM = ^ (p 2 -4M 2 )V( P ,M), (31) 



D 2 (p) = -t^j- 8 {(p 2 + QM 2 )(p 2 -4M 2 ) + 15M 4 }J F (p,M),(m 
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Coefficient 



W 3 
W A 

W 6 



Cut-off regularized value 



N (A 4 -4A 2 M 2 -2M 4 -6M 4 ln (^") 
96 (4tt) 2 /» 
((A 2 +M 2 ) 2 +3M 4 ln(^)) 



N 



(' 



192 (4tt) 2 /» 



-AT ( 15A 2 +2M 2 +30M 2 ln(^j-) 

1440 (4tt) 2 7 s 
-A"(5A 2 -M 2 ) 
960 (4?r) 2 / 8 
AT (l7-60 ln(^)) 
28800 (4tt) 2 / 8 

n (i7-m M*%t)) 

57600 (4tt) 2 / 8 ~ 



Table 1: Regularized constants computed by using a cut-off A. 



with 



J F (p,M) = { 



i-^+i 



(4vr) 2 



(4tt) 2 



2 + 2 ,/4Mi-ltan- 1 | . L 

WES 



27T 



; p 2 <o, 

; < p 2 < AM 2 , 
; AM 2 < p 2 . 



Thus for the particular case of massless branons we find: 



D 1 (p)=2D 2 (p) 



-iNp 4 J F (p, M = 0) Np A ln(-p 2 



480/ 8 



480(4tt) 2 / s 



(34) 



which is in agreement with previous results [9, 10]. 

In the general case, the local actions r^[0] and T^[0] are divergent 
and therefore need to be regularized. By using a cut-off A, the C\ constant 
appearing in rW[0] is given by: 



Ci(A,/) 



iVA 4 
16(4tt) 2 / 4 



(35) 
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The Wi constants can be found in Table 1. A could represent the width of 
brane or any other mechanism that modified the short- distance theory to cure 
the ultraviolet behavior of branons. However, for our purposes, A is just a 
phenomenological parameter. From the point of view of the effective theory, 
A// parameterizes how strongly (or weakly) coupled quantum brane is, and 
therefore controls the unknown relative importance of tree-level versus loop 
branon effects. From (35) we can see that the perturbative loop analysis only 
makes sense for approximately A< A^nf N~ l / A . 

The first term of the effective action T^ 1 ^] is proportional to the trace 
of the energy momentum tensor so it would vanish if the SM were a scale 
invariant theory. By using the equation of motion O^Tg^ = it is possible 
to show that the only effect of this term is to renormalize the boson (scalars 
or gauge fields) masses m$ A = (1 — 2Ci) l ^ 2 mq >i A and the fermion masses as 
m^p = (1 — Ci)m^. Therefore, this first correction to the SM action does not 

have any measurable effect. On the other hand, the local action r^ 2 ' ) [0] is 
important from the phenomenological point of view. At low enough energies 
the dominant terms are the ones proportional to W\ and W^. Thus the SM 
Lagrangian is complemented by the additional effective Lagrangian given by: 

AC eff = W 1 T IX1/ T^ + W 2 T»T». (36) 

where for A > M 

NA 4 



W x = 



96(4tt) 2 / 8 



TV A 4 

w '< = mWT^ (37) 

When this is not the case, one should use the full result in Table 1. From 
this Lagrangian it is possible to obtain the corresponding Feynman rules (see 
Appendix B). The most relevant contributions of branon loops to the SM 
particle phenomenology are the four-fermion interactions and fermion pair 
annihilation into two gauge bosons, whose cross-sections can be found in 
Appendix C. 



4 Phenomenological consequences: constraints 

An effective Lagrangian similar to (36) was obtained in [10, 11] by integrating 
at the tree level the Kaluza-Klein modes of gravitons propagating in the bulk 
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Experiment 


(TeV) 


C (pb- 1 ) 




HERA C 


0.3 


117 


52 


Tevatron-I a ' b 


1.8 


127 


69 


LEP-II a 


0.2 


700 


59 


LEP-II fc 


0.2 


700 


75 


combined 




81 



Table 2: Lower limits on f 2 / (iV^A) (in GeV) from virtual branon searches 
at colliders (results at the 95 % c.L): HERA [12], LEP-II [13] and Tevatron-I 
[14]. The indices a ' b,c denote the two-photon, e + e~ and e + p (e~p) channels 
respectively, s/s is the center of mass energy of the total process, and C is 
the total integrated luminosity. 



and some of its phenomenological consequences were studied there. Thus it 
is easy to translate some of the results from these references to the present 
context. 

Concerning the four fermion interactions ^ a (Pi)'4 ) a(P2) — ► ^bip^ji^bipi) 
(see Apendix B.I.), the most interesting case is the Bhabha scattering at LEP. 
From the Lagrangian in (36), it is possible to find a four- fermion interaction, 
whose amplitude is given by: 

M Alp = v a/ (p 1 )u b (p 4 )V^ w (-p 1 ,-p2,P3,P4)v b/ (p 3 )u a (p2). (38) 

Neglecting the fermion masses and assuming that all of them are different, 
the amplitude is just: 



Wi _ 

M * P (pi,P2,P3,P4) = —Va'{pi)u b {pi) 



la'a^lw(P2 ~ Pl)v(P4 ~ Psf 



+ la'a^lbb'u{P2 ~ P\Y {Pi ~ P<iY 



V b '(p3)u a (p2) ■ (39) 



which can be compared with the analogous amplitude given in literature for 
the graviton case [15, 16]. On the other hand, the explicit form of the cross 
section is given in Appendix C, which agrees with [16, 17]. These interactions 
would lead to modifications of the scattering cross sections such as 

e + e~ -> ££, qq (40) 
qq - U, qq. (41) 

where £ = e, //, r. The Bhaba scattering, fermion pair production in e + e~ 
colliders and Drell-Yan production at hadron colliders have been studied in 
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detail from the point of view of the KK- graviton virtual exchange [16]. Dilep- 
ton and dijet channels have been studied at LEP and dielectron production 
at Tevatron. Also these processes are interesting for e ± p — >■ e ± p interactions 
observed at HERA . 

For V'(pi), i^ipi) — > A a ^{pz), A^ip^) interactions (see Appendix B.2.), 
diphoton production have been studied at LEP and Tevatron, whereas WW 
and ZZ production have been studied also at LEP. Exchange of virtual bra- 
nons can also contribute to processes like 

e+e - , qq -> 77, W + W~ , ZZ and gg (42) 
77, 99 -»• qq- (43) 

The contribution to the gauge boson production is given by 

M 2 ^ 2A = v a ,(pi)V*ilt(-Pi, -P2,P3,P4)u a (p 2 )e^(p 3 , avyfa, a b ) . (44) 

and to fermion-antifermion production: 

M 2A ^ = uAPi)V a 2 ,tl A (Pi,P2, -ps, -Pi)v a (p 2 y(p 3 , v v y(pt, a b ) . (45) 

In fact, for the graviton case, the e + e~ — > / -f-f,W + W~ , ZZ processes 
have been studied in detail in [18], as well as the gg — > Z + /~[16]. Moreover, 
explicit expresions for the fermion-antifermion and diphoton production cross 
sections are given in Appendix C, or alternatively in [17]. 

Using the processes described above, we have obtained limits on the pa- 
rameter combination / 2 /(A^ 1//4 A) from different experiments at LEP, HERA 
and Tevatron, which are summarized in Table 2. Also, with the same anal- 
ogy with the Kaluza-Klein gravitons, we can estimate the constraints from 
future colliders. With that purpose, we have taken into account the estima- 
tions calculated by Hewett [16] for future linear colliders, Tevatron and LHC 
(see Table 3). 

5 Two loops effects: Elect roweak precision 
observables and muon anomalous magnetic 
moment 

Electroweak precision measurements are very useful to find constraints on 
new physics models. The so called oblique corrections (those corresponding 
to the W, Z and 7 two point functions) use to be described in terms of the 
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Experiment 


(TeV) 


C (pb- 1 ) / 2 /(iV 1 / 4 A) (GeV) 


Tevatron-II a ' 6 


2.0 


2 x 10 3 


83 


Tevatron-II a ' fe 


2.0 


3 x 10 4 


108 


ILC 6 


0.5 


5 x 10 5 


261 


ILC 6 


1.0 


2 x 10 5 


421 


LHC 6 


14 


1 x 10 4 


332 


LHC 6 


14 


1 x 10 5 


383 



Table 3: Estimated sensitivities on the parameter / 2 /(A^ 1//4 A) (results in 
GeV at the 95 % c.l.) for various colliders with center-of-mass energies 
and integrated luminosities as indicated. The indices a,b ' c denote the same 
channels as in Table 2. 



S,T,U [19] or the ei,e 2 and e 3 parameters [20]. The first order correction 
coming from the Kaluza-Klein gravitons in the ADD models for rigid branes 
to the parameter: e = dM^/M^ — bM\jM\ was computed in [21]. This 
result can be written as: 

s ^^Mk_m (iWl+2W2) , (46) 

Translating this result to our context as in the previous section we find: 

HMl-M^ N^ 

12(4vr) 4 / 8 K ' 

Notice that this is in fact a two-loop result since it is obtained from a 
one-loop computation by using an effective Lagrangian which is coming from 
another one-loop computation. 

The experimental value of e obtained from LEP [23] is e = (1.27±0.16) x 
10~ 2 . This value is consistent with the SM prediction for a light higgs m H < 
237 GeV at 95 % c.l. On the other hand, the theoretical uncertainties are one 
order of magnitud smaller [20] and therefore, we can estimate the constraints 
for the branon contribution at 95 % c.l. as \5e\ < 3.2 x 10~ 3 . Thus it is 
possible to set the bound: 

f 4 

> 3.1 GeV (95 % c.l.) (48) 



AT1/2 A 3 



This result has a stronger dependence on A (A 6 ) than the interference cross 
section between the branon and SM interactions (A 4 ). Therefore, the con- 
straints coming from this analysis are complementary to the previous ones. 
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A further constraint to the branon parameters can be obtained from the 
H anomalous magnetic moment. The first branon contribution to this pa- 
rameter can be obtained from a one loop computation with the Lagrangian 
given by (36). The result for the KK graviton tower was first calculated by 
[24] and confirmed by [21] in a different way, and can be written as: 

2m 2 A 2 

6a ^ 3(£^ {UWi ~ 12W2) > (49) 

which for the branon case can be translated into: 

5 m 2 AA 6 . . 

^ * Tl4(#— ■ (50) 

This expression is is qualitatively similar to other g — 2 contributions obtained 
in different analyses of extra-dimension models [22]. The result depends on 
the cut-off A in the same way as the electroweak precision parameters. How- 
ever the experimental situation is slightly different. In a sequence of increas- 
ingly more precise measurements, the 821 Collaboration at the Brookhaven 
Alternating Gradient Syncrotron has reached a fabulous relative precision of 
0.5 parts per million in the determination of = {g lx — 2)/2 [25]. These mea- 
surements provide a stringent test not only of new physics but also of the SM. 
Indeed, the present result is only marginally consistent with the SM. Taking 
into account the e + e~ collisions to calculate the 7r + 7r~ spectral functions, the 
deviation with respect to the SM prediction is at 2.6 standard deviations 
[26]. In particular: 5a^ = a^(exp) - a^(SM) = (23.4 ± 9.1) x 1(T 10 . Using 
Equation (50) we can estimate the preferred parameter region for branons to 
provide the observed difference: 

f 4 

6.0 GeV > f /2 > 2.2 GeV (95 % c.l.) (51) 

We observe that the correction to the muon anomalous magnetic moment is in 
the right direction and that it is possible to avoid the present constraints and 
improve the observed experimental value obtained by the E821 Collaboration. 

If there is really new physics in the muon anomalous magnetic moment, 
and this new physics is due to branon radiative corrections, the phenomenol- 
ogy of these particles should be observed at the LHC and in a possible future 
ILC (see Table 3). In particular, the LHC should observe an important dif- 
ference in the channels: pp — > e + e~ and pp — > 77 with respect to the SM 
prediction. The ILC should observe the most important effect in the pro- 
cess: e + e~ — > e + e~. Moreover, in [27] it was shown that the same parameter 
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region in which branons could explain the muon anomalous magnetic mo- 
ment is also compatible with a cosmological branon relic abundance enough 
to account for the observed dark matter relative density [28] . 

6 Conclusions 

In this work we have studied the phenomenological consequences of branon 
radiative corrections, by calculating the one-loop effective action for SM par- 
ticles, obtained after integrating the branon fields out. We have found new 
interaction vertices of SM particles, in particular new four-fermion inter- 
actions and interactions involving two fermions and two gauge fields. Our 
results, computed with a cutoff regulator, are very similar to those obtained 
by integrating the graviton KK modes at the tree level in ADD models, and, 
accordingly, we have translated the different constraints to the branon case. 
Thus, we have obtained limits for the combination of parameters f 2 / (iV^A) 
from present experiments at LEP and Tevatron, and also for future colliders 
(ILC and LHC). 

We have also considered the branon two-loop effect on electroweak pre- 
cision observables and on the muon anomalous magnetic moment. We have 
evaluated the corresponding corrections and obtained the preferred parame- 
ter range for branons in order to fit the Brookhaven results, and at the same 
time, to be consistent with LEP precision measurements. In Fig. 1 we in- 
clude all those limits and also the parameter region in which the theory can 
be considered as strongly interacting, i.e. (A > A^JixfN~ l l 4: ) and for which 
the loop expansion is no longer valid. 
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Figure 1: Main limits from branon radiative corrections in the / — A plane for a model 
with N = 1. The (red) central area shows the region in which the branons account for the 
muon magnetic moment deficit observed by the E821 Collaboration [25, 26], and at the 
same time, are consistent with present collider experiments (whose main constraint comes 
from the Bhabha scattering at LEP) and electroweak precision observables. Prospects for 
future colliders are also plotted. 



A Divergent integrals 

Definitions of various divergent integrals used in the text: 

J(°)(p) = / <fij-2 _ M 2 + i£)((p + q) 2 _ M 2 + ~ ) ■ ( 52 ) 

J$(p) = J dq {q2 _ M2+ u) ^ U + g)2 _ M 2 + u) ■ ( 53 ) 
= [da W»(P + *)p(P + *)>> (M) 

■^vpoKP) J y ^ 2 _ M2 + + g y _ M2 + ie y 

B Effective Feynman rules 

In this section we give the most important effective Feynman rules, obtained 
by the integration of the branons. We are going to present the fundamental 
new vertices with outgoing momenta from (36). 
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B.l Effective 4-fermion vertex 

One of the most relevant contribution of virtual branons to the phenomenol- 
ogy of the SM particles is the effect on four-fermion interactions. For a 
generic four-fermion process, the branons induce a new effective vertex of 
the form: 




Wi r 

= — 32m^ a 5 a > a m^ b 5 w + 7o'a M 7t,b'(P2 - Pi)v{pa ~ PzY 

(P2 - PlY + IZniTpJa'aJbV n(p4 ~ PsY 
+ {la'a^lbb'u + ^la'aulbb' /J G°2 ~ PlY (Pi ~ PzY 

w 2 



+ ~Y [ 8m V>A'a + 3la<aii(P2 ~ PlY] 

[Sm^Sw + ^wv{pa - PzY] ■ 



(55) 



In the case in which the fermion fields a and b are the same, one has to 
take into account two different effects in order to obtain the vertex from the 
above expression. On one hand, a factor of 2 due to the quadratic term in 
the SM energy- momentum tensor is not present, and on the other hand, the 
symmetrization with respect to the change: {pi,P2} <-> {P3,P4} should be 
performed. Therefore, the general form of the four-fermion vertex is given 
by: 



VSu/(pi,P2,P3,Pi) = V2w{Pl,P2,P3,P4) 



4ipD 



(56) 



Vbbt°a(P3,P4,Pl,P2) ~ Kw(Pl>P2,P3,P4) 



4V>£> 
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B.2 Effective ip, ip, A^, A v vertex 

The exchange of virtual branons can also contribute to processes involving 
both a fermion pair and a gauge field pair As shown in the diagram, branons 
induce a new effective vertex, which, with the same momenta assignment as 
before, takes the form: 






&M) 



Va'tlt(Pl,P2,P3,P4) = W X 



[{P2 - P\Y(PZ\P^) ~ (P4,Ps)(P2 - Pl) M ]7a'a v 



+ [(P2 - P\) X {P±\P-5v) ~ (P4,P3)(P2 - Pl)vYla'a n 

~ [(P2 - Pl) X la'a(P4\P3a ~ P3XPia ~ VXa (P4, Ps) )] V 

- [(P2 - Pl)\P^P3u - (P2 - Pl)nPiXP3u 

~(P2 ~ Pl)uP3XP4Ma ( 57 ) 



+ (P2 - Pi) ^1 a' a v + (P2 Pi) via' a fj\ 

- 2 W 2 m 2 A [8m^5 a > a + 37^, a (p 2 - Pi)a]t?, 



pV 



B.3 Effective B p , B a , A^, A„ vertex 

Another effective interaction produced by virtual branon exchange can also 
contribute to processes involving two different gauge field pairs. Indeed, 
branons induce the following effective vertex: 



V™l A ( Pl ,p 2 ,P3,Pl) 



[11/ pa 



4Wi 



P2 fi P3a{PlpPiu ~ PluPip) 



+ PlpP3p(P2aP4v ~ P2uPia) + P3p(P2aPipPlu ~ P2uPiaPlp) 
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B a (p 2 



B p (pi 



AM 



BM 



AM 





AM 



B P (pi) 



MP3 



- (P3, Pa) [P2n(-PluVv P + PlpVua) + Pl^i-pZuVcrp 
+ PlpT]va) + P2uPlpVap + PluP2aVpp) ~ 2PlpP2aVpv] 

~ (Pi , Pa) [P4v (p&rtyip + P3pVpa ) + P3p {.PipVau + Pia Vup ) 

- 2p3pP4uVvp ~ (P4pP3a + KpPfc)v] 

- (P2, Pi) [P3a(PlpVpu - PluV W ) + P3p,PluVpa] 
~ (P2, Ps) [PtolpipTliw ~ VlpVvp) + PAvPlpJlpa] 
~ (P1,Pa) [P3p(P2a11pv ~ P2uV^) + P3pP2uVpa} 
~ (Pi, Ps) [P2 M (P4^p<T - P4p?W) + P4pP2 CT ^] 

- (Pl,P3) ^(jPivVpa - PApVua) + P4pP2aVpv] 
+ (Pl,PA) (P2, Ps) VlpvTlpa + T]ppJ]va - ZVpaVpv] 
+ V P <tVpu[(Pi,Pa)(P2,P3) + (PuP3)(P2,Pa)} 

+ m 2 A [-p 4ll (-2p 3u r] pa + p 3p r] ua + p 3a r] up ) 



~ P3u(pApr] P a + PAaVfip) + ^^(p^pPSc + P3pPAa) 

+ (P3, Pa) (VvaVw + W ~ VpvVpu)] 

+ m 2 B [-p 2p (-2p la r]^ + p lp r] au + PiuVap) 

- Pla {P2pf] pv + P2uVpp) + 2 (j^fiPlv + Pl/iP2v) 

+ (Pi, P2) (^<7^p + ^ VpuVpv)} 



+ m 2 A m 2 B [r] ua r] pp + r]^r] pu ] 



+ 2W 2 m 2 A m 2 B rj pjU rj pu . 



{hi 



In the case in which the gauge bosons A and B are the same, one should 
again take into account the points commented above. In particular the sym- 
metrization with respect to the change: {pi,p 2 , <-> {P3,P4,P, cr} should 
be performed. Therefore, the vertex with four identical gauge boson can be 
written as: 
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Vi„ A pa (PliP2,P3,P4) = 1[V™1 A (P1,P2,P3,P4) + V p lf l/ A ( P3:P4:Pl:P2 )}. (59) 

C Cross sections 

In this section we show the modifications in the cross-section of four-particles 
processes derived from virtual branon exchange in terms of s = (p\ +P2) 2 , 
t = {pi — P3) 2 and u = (j>2 — P3) 2 ■ We are neglecting the masses of these 
particles, which means s + t + u = 0. 

C.l 0-1 : /OO/O2) -> 7b3)7b4) 

The diphoton production by fermion-antifermion annihilation with electric 
charge Qf and number of colours Nf is given by: 



do\ s 2 — 2tu 
~dT ~ N f s 2 



2naQ 2 f W x 



-1 2 



+ — iv 7 ^ 

tu 2 



(60) 



C.2 CT 2 : g(pi)g(p2) -> / + (l?3)/ (P4) 

On the contrary, the dilepton production by gluon annihilation does not 
present interference term: 



da 2 



Wftu 

647TS 2 



2t« 



(61) 



C.3 0-3 : g(pi)g(p2) -> 7(^3)7(^4) 

The situation is similar for the diphoton production by gluon annihilation, 
since there is no SM contribution at tree-level: 



das 
~dt 



647TS 2 



s 4 - 2tu(2s 2 - tu) 



(62) 
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C.4 (74 : e Oi)e+(> 2 ) -> f{pz)f{pi) (f ±v e ,e ) 

To illustrate the four-fermion interaction contribution, we can write the cross- 
section for the fermion-antifermion production (except v e and e~) in e + e _ 
collisions in terms of the vector Vf = Tf — 2Qfsin 2 9 w and axial a/ = Tf 
couplings of the particular fermion field: 



da 4 
~dt 



d<74 

~dt 



+ 



SM 



NfW 2 

1287TS 2 



s 4 - 2tu{bs 2 - mu) 



(63) 



s-Mf 



e f/(t — -u) 3 + a e ajs(s 2 — 6tu) | . 



C.5 cr 5 : e (pi)e + (p 2 ) ->• e b3)e + (^4) 

For the Bhabha scattering, the cross-section presents more terms since one 
has to take into account the t-channel contributions: 



da 5 
dt 



da 4 W[ 

~dt + 1287TS 2 



aWi(Q 
4s 3 I t 



+ 



40s 4 + 6t(31st 2 - 21s 6 - A0s 2 u) + 9t 4 

9s 4 + 22ts 3 + 24t 2 s 2 - llt 3 u - 10t 4 

1 



(64) 



v 2 + a 2 



sin 2 26 w s-M 2 z 



u(4t 2 - As 2 + 5tu) 



+ 



sin 2 26 w t - M\ 



v 2 e (8s 3 + 6ts 2 ) + (v 2 + a 2 )[s 3 + 12s 2 t - 5t 2 (3u + 2t)\ 



All these results are in agreement with the expressions calculated for other 
kind of models which predict the same Lagrangian (36). In particular, for 
KK-gravitons in the ADD model: W 2 = -Wi/(iV + 2) and Wi is related to a 
new energy scale: Ms [16, 17, 18]. For example, W\ = 4A/MJ in [16] (where 
typically A = ±1, takes into account the unknownness of the exact theory). 
So we can use the effective vertices and cross-sections given in the Apendices 
B and C for branons or gravitons using the corresponding definitions of the 
parameters W\ and W 2 - In fact, we can estimate directly the bounds over 
/ 2 /(AA^ 1//4 ) using the bounds over M s . In the most interesting cases, the 
contribution of the term proportional to W 2 is zero or neglegible and it is a 
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good estimation to take, for A = 1: 



f 2 M<? 

- 0.064 M s . (65) 



AATi/4 4(24tt 2 ) 1 /4 

Typically, when / <C Md the most important signal of brane worlds comes 
from branons. In such a case, we can estimate f 2 / (AiV 1 / 4 ) as it is shown in 
Table 2. 
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